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Most conventional techniques for the determination of microalgae Upid content are time consuming and in 
most cases are indirect and require excessive sample preparations. This work presents a new technique that 
utilizes radio frequency (RF) for rapid lipid quantification, without the need for sample preparation. Tests 
showed that a shift in the resonance frequency of a RF open-ended coaxial resonator and a gradual increase 
in its resonance magnitude may occur as the lipids content of microalgae cells increases. These response 
parameters can be then calibrated against actual cellular lipid contents and used for rapid determination of 
the cellular lipids. The average duration of lipid quantification using the proposed technique was of about 1 
minute, which is significantly less than all other conventional techniques, and was achieved without the need 
for any time consuming treatment steps. 

Biodiesel produced from microalgal lipids is considered a promising alternative to fossil fueP For an 
economical biodiesel production from microalgae, high biomass productivity and lipid content are 
required. These two factors are difficult to achieve simultaneously, as conditions favoring high biomass 
productivity usually result in low lipid accumulation, and vice versa^"^. Under stress conditions, wherein growth 
rate drops, microalgae may accumulate larger lipid content^'^. This lipid accumulation should be monitored 
effectively for economical microalgae-biodiesel production. Therefore, a reliable and rapid quantification of lipid 
content is required. Although several conventional methods have been developed to determine the lipid content 
of microalgae ce\W~^^, they are in general invasive, indirect, require time treatment steps and sometimes experi- 
mentally unrepeatable. Depending on the processing steps required, it usually takes from 2 hours to 2 days to 
determine the amount of lipid content^\ On top of that, these methods are inherently not real-time, therefore, 
would be very difficult to use in controlling and rapidly adjusting the growing conditions to achieve a desired 
objective. Therefore, it is desired to have a rapid feedback based monitoring system, which could also be 
integrated in-situ in a control loop, to adjust the growing conditions without the need for any time consuming 
and expensive sample treatment. 

Electrical characterizations using mainly radio frequency (RF) techniques have received increasing attentions 
in the last decades due to their noninvasive characterization ability of unlimited types of materials without the 
need for any sample treatments'^. It has been reported that RF fields that penetrate through microalgae cells are 
affected by the type and composition of the cells'l The RF fields propagated through microalgae suspension are 
altered in magnitude and phase depending on the cells lipid contents and cells concentration. This paper focuses 
on microalgae lipid content determination using RF, which is based on the cells intrinsic material properties that 
can be modeled as smart dielectric materials. Such a technique provides an accurate lipid determination in about 1 
minute, which is significantly faster than the fastest lipid determining technique. A relationship is developed to 
correlate between the scattering parameters (Sn) and the actual lipids content. This relationship was validated by 
comparing its predictions against the actual values. 

Results 

The presented method is based on measuring high frequency scattering parameters (Sn) of microalgae suspen- 
sions placed in an open-ended coaxial RF resonator structure. The microalgae suspension affects the field 
distribution inside the resonator and consequently the input impedance and resonance frequencies are changed. 
The amount of the change depends on the intrinsic physical properties of the microalgae, namely its dielectric 
constant and electrical resistive losses. Since the size of the microalgae used in this work, namely Scenedesmus sp., 
is small, approximately in the range of several of micrometers (Fig. 1(a)), the propagation of RF field polarizes it 
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Figure 1 | (a) Optical photo of the Scenedesmus sp. cell, (b) schematic of the polarized microalgae inside an alternating current electric field. The 

frequency measurements are expected to deviate in magnitude and phase as the lipid content changes. 



due to the induced charging effect (Fig. l(b))^^"^^. The propagated RF 
fields are altered in magnitude and phase depending on the micro- 
algae lipid contents and cells concentration. This is because as the 
signal frequency increases, the charge distribution in the microalgae 
is modulated, which in turn can create a polarizing effect, thus 
enables the lipid quantifications. At a constant concentration of cells, 
the polarization depends on the compositions and size. The micro- 
algae could be presented as an electrical dipole, which has two pair of 
electrical charges of equal magnitude but opposite sign, separated by 
some distance. 

The process of developing a correlation for lipid content as func- 
tion of scattering parameter Si i for a specific kind of microalgae, with 
predetermined cells concentration, starts with taking a sample from a 
bioreactor, diluting it to a constant cells concentration and then 
placing it in a RF resonator structure. A network analyzer is used 
to measure the scattering parameter of the suspended material which 
is recorded and then processed to extract the change in frequency 
and magnitude. The actual lipid contents of the microalgae were 
determined using Soxhlet extraction technique with n-hexane^^ 
and calibrated against their respective RF signature. The n-hexane 
can extract the non-polar lipids only, which are the interesting lipids 
for biodiesel production. Hence, the measurements were limited to 
them only, and the composition is shown in Table 1. 

Details of the processing steps:. Sample preparation. Samples were 
collected at regular intervals and diluted in growth medium to obtain 
a 4 ml of cell suspension. Figure 2(a) shows the bioreactor from 
which the samples were collected over days of starvation. The Nile 
Red stock solution was prepared as described by Siaut et aP^. 
Figure 2(b) shows a comparison between the actual lipid content 
and the fluorescence intensity determined by the plate reader. The 
primary stress applied was nitrogen deficiency, where accumulations 
of more than 70% (dry basis) have been achieved. This is mainly due 
to the lack of nitrogen required for protein synthesis, and the excess 
carbon from photosynthesis is then diverted into lipid production 
pathway. Figure 2(c) shows the fluorescence microscopy images of 
cells stained with Nile Red that shows the lipid accumulation with 
days of nitrogen starvation. 



Table 1 Scenedesmus sp. microalgae 


lipid composition 


Fatty acid methyl ester 


% Composition 


Methyl myristate 


12.5 


Methyl sterate 


6.3 


Cis-9-oleic Methyl ester 


12.5 


Trans -9-elaidic Methyl ester 


43.8 


Methyl linoleate 


18.8 


Methyl behanate 


6.3 



RF Characterization. In order to confirm the concept, a sample of 
Scenedesmus sp. microalgal suspension was withdrawn from the 
bioreactor and concentrated at fixed cells concentration of 1.5 X 
10^ cells/ml, and then inserted inside the open-ended coaxial res- 
onator as shown in Fig. 3(a). The setup was calibrated before each 
measurement has been carried out. This calibration ensured that the 
measurement actually represented the sample under test while 
excluding the effect of losses and phase shifts due to the cables and 
host structure which could add noise to the measured signal. A blank 
was prepared by harvesting the microalgae cells and its lipid content 
was completely extracted, and then rediluted in growth medium in 
the same concentration of 1.5 X 10^ cells/ml. 

The set of Sn parameter measured curves versus frequency for 
different lipid contents are shown in Fig. 3(b). The RF response of the 
suspended microalgae pronounces the lipid accumulations as the 
magnitude of the Sn changes gradually with the increase in lipid 
content. The Sn curve shows a smooth performance over frequency 
and exhibits resonance with specific resonance frequency, band- 
width and magnitude. All of these parameters vary with the lipid 
content. The resonance frequency of samples free of lipids (i.e., all 
lipids have been extracted) was found to be higher than that of 
samples of higher lipid content, due to the change in the effective 
suspended solution dielectric constant. 

Lipid quantification. The lipid content was calibrated against the 
change in resonance parameters. The resonance can be predicted 
with the Gaussian distribution function. Gaussian function para- 
meters were correlated with the RF measured parameters, namely, 
its bandwidth, magnitude and center frequency. A general Gaussian 
function could be expressed by Eq (1): 



G{f)=yo + Ae- 



(1) 



where,/, ^ and/o are frequency, offset, magnitude, bandwidth 
and center frequency of the resonance respectively. These set of 
parameters were correlated against lipid content. The predicted 
and measured Sn response corresponding to lipid contents of 
12.10, 44.5 and 73.2% were drawn against each other as shown in 
Fig. 3(c). The prediction was carried out from 10 GHz to 13.5 GHz, 
the resonance window. Figure 3 (a) reveals a good agreement 
between the predictions and measurements, with a standard devi- 
ation of ±0.001. The set of parameters {yo, A, w and/o} were deter- 
mined for different lipid contents by fitting each measured curve in 
Figure 3(b) with Eq (1). The predicted coefficients are presented in 
Table 2, with each row represents specific lipid content. 

The resonance parameters over lipid content variation are 
depicted in Fig. 4, which shows that the predicted set of parameters 
versus lipid content exhibited a linear relationship. As the lipid con- 
tent increased, the off- set and magnitude parameters changed in the 
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Figure 2 | Conventional lipid determination: (a) bioreactor with outer diameter of 10 cm, an inner diameter of 5 cm and a height of 40 cm. (b) 
comparison between ( O ) lipid content and (•) relative fluorescence intensity, (c) Fluorescence microscopy images of Scenedesmus sp. cells stained with 
Nile Red and showing lipid accumulation after 1, 11, 14, 20 and 23 days of nitrogen starvation in cultures. 



same manner, while both center frequency and bandwidth exhibited 
identical behaviour. 

The model was validated by measuring RF response of a sample of 
unknown lipid content. The measured RF signal was fitted to Eq (1) 
and the associated coefficients were determined and located on the 
calibration curves shown in Fig. 4. The lipid content that correspond 
to these values was then determined to be 12.05%. The actual lipid 
content of the sample was measured using the conventional soxhelt 
technique and found to be of 12.08%, which is almost identical to that 
predicted by the model with an error of only 0.2%. This clearly proves 
that the proposed technique can be used predict accurately the lipid 
content. 

The accuracy of the proposed model is approximately 98%, and 
the recorded error associated with the lower and higher contents did 
not exceed 2% in the entire range of measured lipid content. The 
error can be further optimized by involving more measured points 
and by having micro cavity resonator dedicated design. The cal- 
ibration curve can now be used to determine the lipid content of 
an unknown sample from measuring the Si i frequency response and 
Gaussian prediction. 

Discussion 

New methods for lipid detection and quantification have been place 
since 1993 and are still getting further attention. However, few meth- 
ods managed to scale with the volume or contents. While there is 
tremendous expansion in high algae growth and various envir- 
onmental conditions, numerous researchers are still only providing 
a qualitative analysis of the lipids. Moreover, the time of quantitative 
analysis for these lipids ranges from few hours to days. A method of 
lipid quantification with less than a minute response detection time is 
necessary in controlling and rapidly adjusting the growing condi- 
tions to achieve a desired objective of the growth. 

The RF propagation is confined and guided inside the coaxial 
cable; the open-end of this cable forms a capacitance within which 
the microalgae suspension acts as its dielectric material. The initial 
capacitance corresponds to the lipid free content; exhibits its min- 
imum value with low resistive losses (high quality factor). As the lipid 



content increases the effective dielectric constant of the solution 
increases. Figure 5 shows the extracted capacitance and resistance 
values versus the lipid content. Both capacitance and resistance 
showed a linear relationship versus the lipid content. The change 
in the magnitude level of Sn parameter is associated with the change 
in the DC resistance of the suspended solution, whereas the change in 
the frequency shift is associated with the change in capacitance. Such 
a strategy streamlined the interpretation of our results in terms of 
specific RF signatures for these different lipids contents, which could 
then be attributed solely to the specific nature of respective micro- 
algae particles intrinsic properties. 

The selectivity and the sensitivity are two important issues that 
should be addressed in any detection method. The conventional 
method, used for lipid determination, involves solvent extraction 
that requires separation of the solvent from crude extract after the 
extraction, which is followed by lipid quantification which may 
require fatty acid composition using chromatography tech- 
niques^^"^^. Although these methods have been used widely and have 
a high accuracy in lipid detection, they are time- and labor-intensive, 
making it difficult to screen large numbers of algae. To ensure the 
ability of the detection selectivity of the proposed approach, the 
frequency responses of other strains of microalgae with identical 
concentration have been measured and depicted in Fig. 5(b). 
Figure 5(b) reveals that each microalgae strain was associated with 
a certain shift in resonance frequency and magnitude; i.e. unique RF 
signature. Figure 5(b) shows the RF measurements for two different 
types of microalgae {Scenedesmus sp. And N anno chlorp sis), revealing 
that each strain has its own resonance frequency, which proves 
that the proposed method can be used not only to determine the 
lipid content of other strains, but can also be used to identify the 
strain. The two strains used are different in size, composition, i.e. 
different dielectric constants and resistive losses. The extracted lipid 
shows a higher dielectric constant and better matching as shown in 
Fig. 5(b-iii). 

On the other hand, the sensitivity of the method has been checked 
by diluting a known concentration of Scenedesmus sp. cells with 
cultivation medium solution, and the results are shown in Fig. 3(c). 
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Figure 3 | Measurements setup and results: (a) Radio frequenq^ measurements setup, (b) Sn parameter measurements of suspensions of microalgae 
with different lipid contents and (c) prediction superimposed with measurement. The samples were withdrawn from the photobioreactor shown in 
Fig. 2(a) after different starvation durations. The cells of known lipids content were diluted to a constant cells concentration and then loaded into the 
hosting RF coaxial resonator structure. 



This figure shows the sensitivity of the measurements, which reveals 
that the presented technique is sufficiently sensitive to allow the 
determination of compounds generated during early-stage cell 
growth. The figure shows that the frequency shift decreases as micro- 
algae concentration decreases. The RF screening has a higher select- 
ivity than other conventional methods, and thus allows a better limit 
of detection without the need of any treatment or labeling. 

This work is a basic experimental attempt to demonstrate and 
explore the use of RF method for lipid content determination. The 
conducted measurements and determined calibration model are 
applicable to the specific microalgae strain, grown in the specific 
cultivation medium under certain set of conditions mentioned in 
this work. If the same strain used is grown in the same medium, 
under same conditions, then the effects of any possible relative 



Table 


2 prediction 


coefficients 








days 


Lipid (%) 


yo 




w 


A 


0 


0 


-1.034 


1 1.896 


0.593 


-3.848 


2 


1 1.74 


-0.974 


1 1.897 


0.643 


-3.461 


11 


15.50 


-0.962 


1 1.901 


0.648 


-3.370 


14 


22.93 


-0.938 


1 1.91 1 


0.654 


-3.219 


20 


44.52 


-0.871 


1 1.939 


0.668 


-2.843 


23 


73.27 


-0.775 


1 1 .978 


0.696 


-2.326 



variations are eliminated. The same method can be applied to other 
strains, grown in other media/conditions. To determine the effect of 
all variables affecting the composition of any strain, a complex 
experimental design is required. To confirm that this approach can 
be applied on other strains, a selectivity study has been carried out, as 
shown in Figure 5 (b). 

The proposed method provided a better combination of high 
sensitivity, quick response, of less than a minute, low cost, high 
throughput, and ease of use. When compared with other lipid deter- 
mination techniques, the developed approach was found to be the 
faster than Raman spectroscopy technique, which needs 2 hours, 
and cheaper than NMR. The repeatability studies demonstrate a very 
good accuracy and good predictions. In addition, the proposed tech- 
nique could be upgraded to be applied in-situ, which will not only 
pave the way for direct and rapid determination of microalgae lipid 
content, but also enable the optimization and selection of species and 
growth conditions. Although in this study, this approach was focused 
on the determination of the microalgal lipids content, the analytical 
capability and quantitation algorithms could be applicable to other 
organisms. 

Methods 

Materials. ?i-Hexane, acetone, Nile Red (9-diethylamino-5-benzo[a] 
phenoxazinone), dimethyl sulfoxide (DMSO), were purchased from Sigma- Aldrich 
Inc. 
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Figure 4 | Extracted parameters versus lipid content: (a) the offset, (b) the bandwidth, (c) the center frequency and (d) the resonance magnitude. The 

measured values used for fitting are depicted by ( O ) and the measured point depicted by ( • ) is used for validation. 



Strain and culture medium. Scenedesmus sp. culture was obtained from Algal Oil 
Company, Philippines, and cultivated in a nitrogen deficient modified Bassel medium 
(-N-BBM), which is composed of (in mM); 0.17 CaCl2-2H20, 0.3 MgS04-7H20, 1.29 
KH2PO4, 0.43 K2HPO4, 0.43 NaCl, 1 ml.L"^ of Vitamine B12, and 6 ml.L"^ of P-IV 
solution that consisted of 2 Na2EDTA-2H20, 0.36 FeCl3-6H20, 0.21 MnCl2-4H20, 
0.37 ZnCl2, 0.0084 CoCl2-6H20 and 0.017 Na2Mo04-2H20. 

Microalgae cultivation. The microalgae was grown in a nitrogen deficient medium 
(-N-BBM) for three weeks, to enhance the biomass productivity, in a 5 L bubble 



column photobioreactor with an internal illumination. All cultivations in this work 
were autotrophic, with CO2 naturally present in air bubbled through the system being 
the sole carbon source. 

Prepared medium was sterilized in an autoclave (Hirayama HV-50, Japan) at 
121°C for 15 min and cooled to room temperature prior to use. The photobioreactor 
was illuminated with one 50 cm, 60 watts, white fluorescent light at a light intensity of 
120 |imol.m~^.s"\ measured using a light meter under 12 h light/dark photoperiod 
automatically controlled by 24 hrs timer (S2402, China). The photobioreactor had an 
outer diameter of 10 cm, an inner diameter of 5 cm and a height of 40 cm. 
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Figure 5 | (a) Extracted capacitance (•) and resistance ( O) versus lipid content ratio, (b) selectivity measurements for different microalgae strains (i) 
Scenedesmus sp. and (ii) Nannochlorpsis and (iii) the extracted lipid only, (c) Sensitivity measurements. 
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The cell growth was monitored daily by measuring the optical density at 680 nm 
using a spectrophotometer (Shimadzu UV-1800 UV, Kyoto, Japan). The samples 
were measured twice and the average values were considered. The cells concentration 
(cells.ml"^) at any given cultivation time was calculated from a pre-prepared cal- 
ibration curve of the optical density at 680 nm versus cells concentration determined 
using Neubauer Hemocytometer, placed on a microscope (Eclipse LVlOO Pol, Nikon, 
Japan). The dry weight of algal biomass was also determined by filtering the algal 
suspension using a Whatman filter paper, dried overnight at 105°C in an oven 
(Memmert, Germany) until constant weight. 

Calibration curve procedure. Cells concentration was determined by the cell 
counting technique, where 0.1 mm deep hemocytometer (PZO, Polandcountry). 
Different culture dilutions were prepared to calibrate cells concentration versus 
adsorption using UV spectrophotometer (Hach DR-5000, Germany) at 680 nm. Two 
ml of the well mixed culture withdrawn from the flask and diluted to 4 ml in growth 
medium. This dilution step was repeated four times. 0.1 ml droplet of 1 : 16 diluted 
were placed on a slide of 0.0025 m^, using Pasteur pipette, and covered with 22 X 
22 mm square cover slip. Cells were allowed to settle for 2 min and then the number 
of cells was counted in the four squares at the corner of the hemocytometer grid of 
nine large squares of 1 mm^ surface area. Total cells concentration per ml was 
calculated using Eq (2), from which the initial concentration was determined by 
multiplying the concentration by 16. 



Cells concentration = 



Total counted cells 



xlO'^ 



(2) 



A calibration curve between cell count and the spectrophotometer absorption reading 
was generated and used in subsequent analysis to convert the spectrophotometer 
readings into cells concentration. The cells counting and dry weight analysis are 
standard procedures to determine cells and biomass consecrations. Further details of 
the procedure can be found elsewhere^°'^\ 

Lipid content determination. Samples were collected at regular intervals and diluted 
in growth medium to obtain a 4 ml of cell suspension. The accumulations were 
monitoredby staining the constant concentration samples (1.5 X 10^ cells.ml"^) with 
Nile Red that emits a yellow fluorescent signal in the presence of the lipid, and the 
fiuorescents were visualized using fluorescence microscope (Olympus). The Nile Red 
stock solution was prepared as described by Siaut et al" by dissolving 0.1 mg of Nile 
Red in 1 ml acetone, and the solution was stored in the dark at 4°C. Culture samples 
(500 |il) were placed in an eppendorf tube, span in a centrifuge (Sigma 113, Germany) 
for 30 s at 4,000 rpm and 410 |il of the supernatant were taken. DMSO (10 was 
then added to promote the accessibility of Nile Red into the cells. The culture was then 
vortexed and 1 |il of Nile Red solution was added followed by 20 min incubation in 
the dark. The lipid accumulations were then quantified using Multi-label Plate Reader 
(Perldn-Elmer, Boston) with black 96-well plates. The cells were harvested by 
centrifugation and then lyophilized in a freeze drier (Telstar, Terrassa, Spain) 
operated at — 54°C and 0.02 mbar for 6 hrs. The actual lipid content was then 
determined using n-hexane in Soxhlet system. 

Radio frequency measurements. The network analyzer (NA) from Rohde and 
Schwartz (R&S®ZVL) is commonly used to measure S-parameters, because reflection 
and transmission of electrical networks are easily measured at high frequencies. The 
NA involves radio frequency generators that produce power signals at different 
frequencies, thus when these power signals propagate inside the cultivation medium 
they scatter back towards the source and some portion of this signal power is absorbed 
by the microalgae particles. The system was calibrated using the SLOT transmission 
line techniques for the network analyzer. A typical calibration moves the 
measurement reference planes to the end of the test cables. Therefore, it excludes the 
effect of losses and phase shifts that could add noise to the measured signal. Each 
microalgae suspension, of different lipid content, was exposed to a radio frequency 
signal with a power of 10 dBm and with a sweep from 8 GHz up to 13.6 GHz (the 
equipment measurements range capability). The self- resonance frequency of the 
coaxial cables was ensured to be pushed above 30 GHz; so it does not affect 
the measurements in the mentioned range. 
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